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We study theoretically quantum melting transitions of stripe order in a metallic environment, 
and the associated reconstruction of the electronic Fermi surface. We show that such quantum 
phase transitions can be continuous in situations where the stripe melting occurs by proliferating 
pairs of dislocations in the stripe order parameter without proliferating single dislocations. We 
develop an intuitive picture of such phases as "Stripe Loop Metals" where the fluctuating stripes 
form closed loops of arbitrary size at long distances. We obtain a controlled critical theory of a 
few different continuous quantum melting transitions of stripes in metals . At such a (deconfined) 
critical point the fluctuations of the stripe order parameter are strongly coupled, yet tractable. They 
also decouple dynamically from the Fermi-surface. We calculate many universal properties of these 
quantum critical points. In particular we find that the full Fermi-surface and the associated Landau 
quasiparticles remain sharply defined at the critical point. We discuss the phenomenon of Fermi 
surface reconstruction across this transition and the effect of quantum critical stripe fluctuations on 
the superconducting instability. We study possible relevance of our results to several phenomena in 
the cuprates. 



I. INTRODUCTION 

In the last several years evidence for the occurrence of 
stripe and related orders has accumulated in many under- 
doped cuprates. Static charge and spin stripes have long 
been known to occur in some La-based cuprates^. Other 
cuprates have been thought to have incipient stripe or- 
dering ( "dynamic stripes" ) which can be pinned locally 
around impurities or near vortex cores in the low temper- 
ature superconducting state. Such pinning of incipient 
stripe order is routinely observed in STM experiments 
in the Bi-based^ and oxychloride 3 cuprates. Further a 
number of phenomena in the YBCO family have indi- 
cated the possible role of incipient stripe order. Quantum 
oscillations^ seen in high magnetic field and low tempera- 
ture have been interpreted^ in terms of the reconstruction 
of a large Fermi surface by stripe order. Such a Fermi 
surface reconstruction might also explain the sign of the 
Hall effect 6 and various other transport properties of un- 
derdoped YBCO at low temperature and in high mag- 
netic fields^. Several experiments also show that YBCO 
undergoes a rapid crossover to a regime of enhanced or- 
thorhombicity at a temperature that roughly tracks the 
pseudogap temperature^—. This is reasonably associ- 
ated with enhanced electronic nematic order which is nat- 
urally a precursor of stripe order. Most recently direct 
evidence for charge stripes in YBCO has been obtained 
in an NMR experiment in high magnetic fields^. 

The stripe ordering seems to disappear with increas- 
ing doping. The Fermi surface of overdoped Tl-2201 has 
been mapped out in great detail through angle dependent 
magnetoresistance studies^, quantum oscillations^, and 
angle resolved photoemission experiments^. All these 
probes clearly and convincingly demonstrate the exis- 
tence of a large band structure-like Fermi surface at low 
temperature. The absence of any Fermi surface recon- 
struction strongly suggests the absence of stripe order- 
ing. 



These developments sketch a picture of the evolution 
of the low temperature "underlying normal" state elec- 
tronic properties upon moving from the overdoped to 
the underdoped side. Decreasing doping tends to pro- 
duce stripe ordering which reconstructs the large Fermi 
surface. It is natural then to expect the existence of a 
quantum phase transition between the large Fermi sur- 
face metal and a stripe ordered metal with reconstructed 
small Fermi pockets. Such a transition, if second order, 
might conceivably play a role in the mysterious phenom- 
ena characterizing the strange metal regime of optimally 
doped cuprates. This viewpoint is advocated e.g. by 
Taillefer in Ref. GJ. 

Despite this strong motivation, there is very little theo- 
retical understanding of such quantum phase transitions. 
In a weakly interacting Fermi liquid the stripe order may 
be viewed simply as a uni-directional charge density wave 
(CDW) or spin density wave (SDW). The transition to 
this kind of order may then be described in terms of a 
fluctuating order parameter coupled to the particle-hole 
continuum of the metallic Fermi surface. Such a theory 
was formulated by Hert a 16 ' 17 and others in the 1970s and 
has received enormous attention over the years. Despite 
this the theory is very poorly understood in two space di- 
mensions (the case relevant to cuprates) . Very interesting 
recent work shows that the low energy physics involves 
strong coupling between the various gapless degrees of 
freedom - the resulting theory currently has no controlled 
description and remains to be understoodiSr— . 

In light of this, and in light of the fact that the cuprates 
are in any case unlikely to be correctly described as 
weakly interacting Fermi liquids, it is natural to explore 
alternate strong coupling approaches to phase transi- 
tions associated with stripe ordering^ and the associated 
Fermi surface reconstruction. Specifically it is interest- 
ing to view this transition as a melting of stripe order by 
quantum fluctuations. Could the stripes melt through a 
continuous quantum phase transition? What is the na- 
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ture of the resultant melted phase? How does the metallic 
Fermi surface affect (and is affected by) such a putative 
continuous stripe melting transition? What is the correct 
description of the universal singularities associated with 
such a stripe melting quantum critical point? 

In a recent short paper we initiated a study of these 
questions with the main goal of explaining some old neu- 
tron scattering experiments in the cuprates. In this paper 
we expand in detail the ideas on stripe melting we dis- 
cussed in our earlier work-22. We provide a few concrete 
and tractable examples of continuous stripe melting tran- 
sitions in the presence of a metallic Fermi surface. As ex- 
pected the stripe melting is accompanied by a reconstruc- 
tion of the metallic Fermi surface. Remarkably despite 
this right at the quantum critical point the critical stripe 
fluctuations decouple dynamically from the particle-hole 
continuum of the Fermi surface and are described by a 
strongly coupled though tractable quantum field theory. 
In the language of renormalization group theory the cou- 
pling of the stripe fluctuations to the Fermi surface is a 
"dangerously irrelevant" perturbation - though it is ir- 
relevant at the critical fixed point it is relevant at the 
stripe ordered fixed point and leads to the Fermi surface 
reconstruction. One consequence of this dangerous irrele- 
vance is that close to the quantum melting transition, the 
energy scale associated with the onset of stripe order is 
parametrically larger than the energy scale at which the 
Fermi surface reconstructs. We determine the universal 
critical singularities both of the stripe order parameter 
and of the electronic excitations at the hot spots on the 
Fermi surface that are connected to each other by the 
stripe ordering wavevector. 

The stripe melting transitions discussed in this pa- 
per and in our earlier work-2 2 . are obtained by prolifer- 
ating some but not all topological defects of the stripe 
order parameter. The resulting stripe melted phase re- 
tains a memory of the long range stripe ordered phase 
by possessing gapped excitations associated with the un- 
proliferated topological defects. Phases of this kind were 
proposed by Zaanen^ and co-workers and further ex- 
plored in Refs. I2414261 Despite having the same sym- 
metries, these are not regular Fermi liquids, due to the 
unproliferated defects. A clear distinction between these 
non-trivial stripe liquid phases and the regular Fermi 
liquid lies in the topological structure, i.e. the former 
have ground state degeneracies on non-trivial manifolds. 
While this difference is completely sharp theoretically, it 
is very difficult to detect with any conventional experi- 
mental probe. In particular these phases have large Fermi 
surfaces with Landau quasiparticle excitations described 
within the usual Fermi liquid theory paradigm. Their 
single-electron properties, transport and low-energy ther- 
modynamics are Fermi-liquid like, thus such phases may 
have easily been mistaken for Fermi liquids. The topolog- 
ical structure associated with the unproliferated defects 
leads to a fractionalization of the stripe order parameter 
- this is extremely difficult to probe in experiments. 

We provide a simple physical picture of the stripe flue- 
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FIG. 1: Left: In the Stripe Loop Metal phase, only those 
stripe fluctuations that result in closed loop patterns are al- 
lowed. Right: In regular fluctuating stripes, both closed and 
open patterns are possible. 



tuations that lead to this kind of stripe melted phase. 
When the stripes melt the resulting ground state will in 
general be a quantum superposition of arbitrary stripe 
configurations. The fluctuating stripe phases described 
in this paper may be pictured as ones in which the stripes 
form closed loops of arbitrary size while they fluctuate 
(see FIG.Q]). Cutting a stripe loop open to leave an open 
end for a stripe costs a finite energy and describes an 
excited state. Such an open end is precisely the gapped 
unproliferated topological defect that survives the stripe 
melting transition. This picture is justified and elabo- 
rated in detail below. We thus dub such phases "Stripe 
Loop Metals" . In contrast the conventional Fermi liquid 
may be viewed as consisting of fluctuating stripes of all 
kinds including ones with open ends in its ground state. 

In the early work such a Stripe Loop Metal was sug- 
gested as a candidate for the underdoped cuprates. As 
it seems extremely unlikely that the "underlying" nor- 
mal state in the underdoped cuprates has a large un- 
reconstructed Fermi surface the Stripe Loop Fermi liq- 
uid is unlikely to occur in the underdoped side. It is 
more interesting therefore to explore the possibility that 
it may describe the "normal" ground state of the over- 
doped cuprates. Indeed none of the existing experimen- 
tal probes of the overdoped normal state are in a position 
to distinguish between such a stripe-fractionalized Fermi 
liquid and the conventional Fermi liquid. Further as we 
demonstrate in this paper the Stripe Loop Fermi liquid 
admits a direct and interesting second order transition 
to the stripe ordered metal with a reconstructed Fermi 
surface. However our results also demonstrate that the 
corresponding critical stripe fluctuations cannot by them- 
selves account for most of the observed non- Fermi liquid 
physics around optimal doping in the cuprates. Never- 
theless as we discuss, this kind of stripe melting transition 
may contain interesting lessons to at least understand the 
nature of stripe fluctuations near optimal doping. Quan- 
tum melting transitions between the stripe ordered metal 
and the conventional Fermi liquid are more challenging 
to address, and we leave them for the future. 

Apart from developing an intuitive description of the 
Stripe Loop Metal phase, we also discuss several experi- 
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mental consequences of the theory of the quantum phase 
transition to the stripe ordered metal. The availability of 
a theory of a continuous stripe melting transition enables 
us to reliably address the phenomenon of Fermi surface 
reconstruction across such a transition. For instance we 
determine how the gap opens at the hot spot as the tran- 
sition is approached from the stripe melted side and the 
growth of the quasiparticle weight in the folded portion 
of the reconstructed Fermi surface in the stripe ordered 
phase. We also study the issue of the low temperature su- 
perconducting instability of the metal and its interplay 
with the stripe melting. Not surprisingly we find that 
the energy scale of the superconducting instability is en- 
hanced as the stripe quantum critical point is approached 
from either side. 

The rest of the paper is organized as follows: In Section 
HI! we discuss some of the challenges presented to theory 
by well known experimental results^!. We then set the 
stage for discussing phase transitions by first describing 
various possible stripe order parameters that are of rel- 
evance for the cuprates in Section [ITU In Section HVl we 
discuss possible phase diagrams for stripe melting tran- 
sitions and the associated defects in the order parame- 
ters. The concept of "Stripe Loop Metals" is introduced 
in Section [V] and connected to the corresponding quan- 
tum field theories. In Section IVII the theories of several 
charge-stripe melting transitions are presented in detail. 
Spin-stripe melting transitions are discussed in Section 
IVHI In Section I Villi we calculate single particle prop- 
erties close to the phase transitions. The possibility of 
pairing by stripe fluctuations within our theory is ana- 
lyzed in Section ITXl In Section [X] we review some results 
of scaling theory applied to fluctuating stripes. We dis- 
cuss some existing experimental results as well as predic- 
tions of our theory for future experiments. Some more 
technical details can be found in the appendices. 



II. CHALLENGES FROM PRIOR 
EXPERIMENTS 

If quantum criticality associated with onset of stripe 
order is held responsible for the physics of the strange 
metal in the cuprates, then it is very important to exper- 
imentally establish that critical stripe fluctuations occur 
in the strange metal regime. There is actually very little 
information from experiments on critical stripe fluctua- 
tions in near optimal cuprates above their superconduct- 
ing transition termperature. 

In an important and well-known experiment, Aeppli et. 
al£l measured the dynamic spin-susceptibility of slightly 
underdoped (x = 0.14) LSCO near (n, 7r) over a wide 
range of temperatures and wave- vectors. However as em- 
phasized in our previous work22 the results paint a very 
intriguing picture of the stripe fluctuations and pose a 
challenge to theory. From the scaling of the width of the 
neutron scattering peak (the inverse correlation length) 
with temperature they deduced that the dynamical expo- 



nent z ~ 1. They further found that for low frequencies 
the peak height scales as T~ 2 . Within z = 1 scaling, 
a standard scaling argument shows that this implies an 
anomalous exponent r\ rs 1 for the critical spin fluctua- 
tions (see Sec. IX Ap . Such a large value of the anomalous 
dimension r\ 1 of the spin-fluctuations is very unusual 
for a Landau quantum critical point, but it is quite com- 
mon for Non-Landau quantum critical points^— . Later 
on in this paper we provide a careful discussion of the 
data of Ref. l27| and highlight the need for further exper- 
imental studies. 
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FIG. 2: The Fermi surface of LSCO at x — 0.15 as measured 
in Ref. [35| via ARPES, shown in the first Brillouin zone. 
Left: Hot-spots on the Fermi surface that are connected by 
the SDW wave vectors. Right: Hot-spots of the CDW wave- 
vectors. 

We would like to emphasize that in a metal, z = 1 
is very surprising. Since the stripe-ordering vector Q 
connects two points on the Fermi-surface (see FIG. 0, 
the stripe-fluctuations should be Landau damped. The 
propagator of the stripe fluctuations is expected to be 
modified as (see FIG. [3]) 

xstripe (fc,a,)= -^— J^_^ (1) 

Landau Damping 1 

oj 2 + ijui — v 2 (k — Q) 2 ' 

where the damping rate 7 is determined by the proper- 
ties of the conduction electrons close to hot spots, i.e. 
points on the Fermi-surface which are connected by Q. 
For small frequencies lo — > the quadratic term can be 
neglected against the damping term and the stripe fluctu- 
ations become strongly coupled to the Fermi surface. In 
particular, we see that oj ~ 5k 2 , where Sk = (k — Q), 
i.e. the dynamical exponent z = 2. Recent careful 
analyses^— show that at low energies the stripe fluc- 
tuations couple even more strongly to the Fermi-surface, 
and higher order contributions become important. There 
is currently no controlled description of the resulting the- 
ory. 

This form of the Landau damping is not unique to 
the weakly interacting Fermi liquid but rather general 
in the presence of gapless excitations into which a stripe 
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fluctuation can decay. In Appendix [5] we demonstrate 
this explicitly for several known non-Fermi liquid metals. 

The data of Ref. [27] shows dynamical exponent of 
2 = 1 over all measured temperatures, i.e. no evidence 
of the expected strong coupling to the Fermi surface is 
seen. Of course, once a single-particle gap develops, a 
low-energy stripe fluctuation can no longer decay into 
particle-hole excitations and effectively decouples from 
the Fermi-surface. This does not require phase-coherence 
and is therefore already possible in the pseudogap phase. 
While such an opening of a gap would explain the ob- 
served z = 1 below the pseudogap temperature T*, it 
would also predict a dramatic change in the nature of 
the stripe-correlations upon crossing T*. Above T* the 
stripe fluctuations are strongly coupled to and modified 
by the gapless Fermi surface, while below T* they de- 
couple. However in the data of Ref. H3, the measured 
spin-correlations develop smoothly across the pseudogap 
temperature, which for x — 0.15 LSCO is around 150 
K2£. 
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FIG. 3: Landau damping of the stripe fluctuations is given 
by the fermionic polarization function near the ordering wave- 
vector. 

The apparent message from the experiments is that the 
critical stripe fluctuations are indifferent to the fate of 
the electronic Fermi surface. This has also been empha- 
sized by other neutron studies of the cuprates where mag- 
netic scattering near the incommensurate peaks seems to 
not know about the particle/hole continuum of the Fermi 
surface at low energies^!. This state of affairs in exper- 
iments should be contrasted with the emerging picture 
from modern clarifications^— of the standard weak cou- 
pling "Hertz-Millis" approach to stripe ordering quantum 
phase transitions: In this theory the stripe fluctuations 
become more and more strongly coupled to the Fermi 
surface at low energy. 

Given this stark contradiction, we therefore focus our 
attention on an alternate theory of the stripe ordering 
transition which describes it as a quantum melting of 
stripes by proliferation of topological defects. 



III. POSSIBLE STRIPE ORDER PARAMETERS 

Throughout the paper we will use the term "stripe" 
as a collective term for various kinds of charge density 
wave (CDW) and spin density wave (SDW) states. With 
the cuprates in mind, we focus on two-dimensional sys- 
tems (a single copper-oxygen plane) with an underlying 



orthorhombic or tetragonal lattice. The simplest stripe 
pattern is a 

1. unidirectional charge stripe. Here the expectation 
value of the charge is modulated at some wavevec- 
tor Q c (see FIG. EJ}, i.e. 

P(r) = (4,aCr,«r) = Po + [Pq/ Q ^ + CC.) , (2) 

where c\. a creates a spin a electron at r and sum- 
mation over spin indices is implied. Such unidirec- 
tional stripe patters occur naturally in orthorhom- 
bic crystals, where there is a preferred direction for 
Q c . In tetragonal crystals unidirectional order is 
possible by spontaneously breaking the lattice ro- 
tation symmetry (in addition to the lattice trans- 
lation symmetry along Q c ). This case is frequently 
referred to as smectic in the literature. 

When the ordering vector Q c is commensurate with 
the lattice, i.e. Q c — — (m x ,m y ), where m IjS are 
integers, this is referred to as a phase with com- 
mensurate charge stripes. For the cuprates, com- 
mensurate stripes with m y (m x ) = and a period 
of m x (m y ) = 4 are most commonly observed ex- 
perimentally. 

2. In a unidirectional spin stripe the expectation value 
of the spin undergoes spatial modulations, i.e. 

S(r) = = e iQ *' r M + e~ iQ >- r M* . (3) 

For collinear spin order, i.e. S(r) x S(r') — 0, it 
follows that 

M = e iB N, (4) 

where A is a real vector. A further possibility is 
spiral order, i.e. 

M = ni + in2 (5) 

where h*\ and are real vectors with ni • n,2 = 0. 
With the cuprates in mind, we consider collinear 
spin-stripes exclusively. In the cases of interest 
here, SDW order at Q s will be accompanied by 
CDW order at Q c = 2Q S . In a weak coupling, 
Landau- Ginzburg approach, this follows since the 
charge density has the same symmetries as the 
square of the spin density. Since the ordering wave- 
vectors of charge and spin are tied together, it 
is sufficient to name either one. We will adopt 
the convention that the period of a commensurate 
stripe refers to the charge sector, i.e. a "period-m 
stripe" has period m for the charge, but period 2m 
for spin (see FIG. I3J. It is often convenient to ex- 
press a spin-configuration not in terms of the spin 
directly, but in terms of the Neel- vector (see FIG. 
H]). In particular, in the experimentally observed 
antiphase stripes, the Neel-vector changes sign be- 
tween two adjacent charge stripes. 
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FIG. 4: A period-4 charge stripe accompanied by an antiphase 
spin stripe. The sign of the Neel-vector changes from one 
charge-stripe to the next. 



On a tetragonal lattice, in addition to the above, 
a checkerboard pattern which respects the lattice 
rotation symmety is possible. Here the charge and 
spin-densities have equal modulations in both the 
x and y directions, i.e. 



(p(r)) = pa + (p Qc 



JQc 



+ c.c. 



(6) 



and likewise for the spin-density. All of these order- 
ing patterns are of relevance in the cuprates, and 
we will discuss them in turn. 



IV. QUANTUM MELTING OF STRIPES: 
POSSIBLE PHASE DIAGRAMS 



Before addressing the quantum phase transition in de- 
tail, we will now briefly discuss a possible phase diagram 
and some properties of the phases which we consider. 
The limiting cases are, on the one hand, the usual Fermi- 
Liquid with a large Fermi surface which respects all sym- 
metries of the underlying crystalline lattice. On the other 
hand, there are various striped phases where translation 
symmetry is broken by static spin and charge order as dis- 
cussed in the previous section. In these phases the orig- 
inal large Fermi-surface is reconstructed. As a striped 
phase undergoes a transition (or a sequence of transi- 
tions) into the Fermi-Liquid, the symmetries broken by 
stripe-order are restored. Depending on the stripe order 
in question, these are spin-rotation, lattice rotation and 
lattice translations symmetry. In general, there may be 
intermediate phases where only a subset of the symme- 
tries are restored. For example, when spin-rotation sym- 
metry is restored but lattice symmetries remain broken, a 
spin-striped phase turns into a phase with charge-stripes 
only. Intermediate phases without stripe-order are also 
possible, e.g. a spin-nematic phase where lattice sym- 
metries are restored but spin rotation symmetry remains 
broken. 



A. Structure of defects in the stripe order 
parameter 

Just like melting of an ordinary crystal, the melting 
of stripes is most conveniently described in terms of de- 
fects in the stripe order. In the limiting case of per- 
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FIG. 5: Unidirectional and Checkerboard stripe patterns. 
The gray lines denote minima in the charge density and the 
arrows denote the orientation of the Neel vector. 



feet long-range stripe order, no defects are present. The 
Fermi liquid lies in the opposite limit, where stripe order 
is completely destroyed and the number of defects is no 
longer well defined, i.e. it may be viewed as a condensate 
of all possible defects. Intermediate phases can be real- 
ized (and characterized) by proliferating only a subset of 
allowed defects. 

Before addressing the phase transition, it is necessary 
to identify the possible topological defects. We begin 
by considering incommensurate, unidirectional spin and 
charge stripes. To allow for defects we must allow the or- 
der parameters to vary in space and time. It is convenient 
to introduce the phases of the stripe-order parameters as 



e s , c (r,t)= [Q s Jr,t)-Q s jO,0)]-r. 



(7) 



9(r, t) has the interpretation of being (Q times) the local 
displacement of the stripes in the Q-dhection at time t. 
The fundamental topological defects are a) dislocations 
where the phase of the spin-order winds by tt and N ro- 
tates by 7T and b) dislocations where the phase of the spin 
order 6 S winds by 2ir (see FIG. [5]). Around these defects 
the phase of the charge stripe 9 C = 26 s winds by 2tt and 
4-7T, respectively. We will hence refer to defects of type a) 
as single dislocations and to defects of type b) as double 
dislocations. Spin order is frustrated around single dislo- 
cations but not around double dislocations (see FIG. IB]). 
In the commensurate case of a period m charge-stripe 
Q c = we can also have elementary oriented do- 

main walls where the entire stripe-pattern is shifted by 
one lattice spacing (9 C increases by ^). A single dislo- 
cation is then located at the point where m such domain 
walls meet (see Fig. EJ. 

In a conventional (weak coupling) approach, the stripe 
melting is described as an ATF-transition for the charge- 
stripe order parameter e' lSc . At this transition, all dis- 
locations proliferate and the resulting state is the Fermi 
liquid. However, in the presence of a Fermi surface, this 
approach becomes problematic for the reasons discussed 
in Section|Hj It is then natural to ask whether these com- 
plications are still present at a modified transition where 
only certain dislocations proliferate. 

A physical mechanism that may lead to a situation 
where some, but not all defects proliferate was already 
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FIG. 6: Left: Single dislocations in the charge stripes are 
bound to half-dislocations for the spin-order, leading to frus- 
tration. Right: Double dislocations in the charge stripes avoid 
frustration. The gray lines denote minima in the charge den- 
sity and the arrows denote the orientation of the Neel vector. 




FIG. 7: Left: Four elementary directed domain walls (dashed 
lines) meet at a stripe dislocation. The charge order in neigh- 
boring domains is related by a primitive lattice translation. 
The thin vertical lines are drawn as a guide for the eye, and 
are separated by one lattice spacing. Right: Configuration 
with three domain walls but without a dislocation. 



mentioned above. We have seen that when spin order is 
present, it is frustrated around single dislocations. This 
raises the energy cost of such dislocations^, potentially 
even above the energy of double dislocations which are 
not frustrated. Thermal melting by proliferating double 
dislocations was studied for spin stripes in Ref. 39 and for 
a striped superconductor/the Larkin-Ovchinnikov (LO) 
state in Refs. l4fj||4ll Quantum melting of the LO state 
was discussed in a cold-atoms context in Ref. |42|. 

Even in the absence of long-range spin order, strong 
short-range spin correlations may significantly raise the 
core-energy of single dislocations. In this paper we thus 
explore a stripe melting transition where double disloca- 
tions proliferate, while single dislocations remain gapped. 
In the resulting stripe liquid phase, all lattice symmetries 
are restored. This phase was first proposed byZaanen^ 
and co-workers and further explored in Refs. |24 - |26| . De- 
spite having the same symmetries, it is not the regu- 
lar Fermi liquid due to the gapped single dislocations. 
A clear distinction between the non-trivial stripe liquid 
phase and the Fermi liquid lies in the topological struc- 
ture, i.e. the former has a ground state degeneracy on 
non-trivial manifolds. While this difference is completely 
sharp theoretically, it is extremely difficult to probe ex- 
perimentally. In particular, the stripe liquid phase has 
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FIG. 8: Schematic zero-temperature phase diagram close to 
the multicritical point X which separates the four phases dis- 
cussed in the text. The dashed line is parametrized by g. 



the same large Fermi surface with well denned Landau 
quasiparticles as the weakly interacting Fermi liquid. Its 
single-electron properties, transport and low-energy ther- 
modynamics are also identical, thus this phase may have 
easily been mistaken for a Fermi liquid. A detailed dis- 
cussion of this phase along with a more physical picture 
is given in Section fVl 

Above we saw that the order parameters b — e l9s and 
N are intertwined only around single dislocations, while 
around double dislocations they are independent. Thus 
at energies well below the core energy of single dislo- 
cations, b and N become separately well denned (this 
will be explained more formally in Section IVBp . We 
can envisage four different phases (see FIG. |5J|, in all of 
which the single dislocations are gapped. First, (N) 7^ 0, 
(b) 7^ is the ordered phase with both spin and charge 
order. Second, (N) = 0, (b) ^ describes a phase with 
a charge stripes but without spin order. Third when 
(N) 0, (b) = 0, lattice symmetries and time reversal 
symmetries are restored. Thus in this phase the spin ex- 
pectation value S = but spin rotation symmetry is bro- 
ken by a spin quadrupole moment Q a 6 = N a Nb — \N 2 8 a b 
(such a phase is frequently referred to as spin nematic). 
Finally (b) = 0, (N) = describes the non-trivial stripe 
liquid we mentioned above and which we discuss in detail 
in the following section. 

In the presence of itinerant fermions the two latter 
phases have a conventional large Fermi surface while in 
the striped phases with (b) =/= the Fermi surface is re- 
constructed by the stripe order. 



V. STRIPE LOOP METALS 

A. Physical picture 

As we mentioned above, the distinction between the 
Stripe Liquid Metal obtained after proliferation of paired 
dislocations and the regular Fermi-Liquid is topological 
and thus difficult to detect. On a pictorial level, how- 
ever, the difference between these phases becomes very 
intuitive, by looking at the charge density ~ e l< ^ a ' r b 2 . 
In a perfectly stripe-ordered phase all the stripes extend 
across the entire system, and there are no stripe endings. 
In the regular Fermi liquid, all kinds of fluctuations are 
allowed, in particular there are dislocations where stripes 
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end (see FIG. [TJ. We are instead interested in a stripe 
liquid phases where only double dislocations are allowed. 
Introducing any number of (static) paired dislocations 
into a stripe ordered phase turns the stripe pattern into 
a set of closed stripe-loops. In the stripe liquid phase, 
where these dislocations proliferate, there will thus be a 
fluctuating pattern of stripe-loops, but no stripe endings 
(see FIG. H}. 

Indeed, this picture can be made precise by connecting 
the stripe patterns to features of the field theory, which 
we will demonstrate below. To see the topological nature 
of the SLM it is useful to consider the system on a cylin- 
der, around which the stripes wind in the ordered phase. 
In the SLM phase, all fluctuations of stripes are allowed, 
as long as there are no loose ends. Thus the number 
of closed charge stripes winding around the cylinder is 
conserved mod 2 (see FIG. [9]). There are two ground 
states which are distinct by the parity P c of closed stripes 
around the cylinder. 




FIG. 9: Three smoothly connected configuration of the Stripe 
Loop Metal on a cylinder. The number of stripes winding 
around the cylinder is conserved mod 2. The parity of closed 
stripes is a topological quantum number distinguishing be- 
tween two possible ground states. 

The distinction between these two states is clearly 
topological, i.e. cannot be determined by just looking 
at any finite region of the cylinder. 



B. Gauge theory 

To make the Stripe Loop representation of these phases 
more precise we now connect it to the corresponding field- 
theory. Formally to express the physical order parame- 
ters in terms of b and N, i.e. 

M = Nb (8) 
ip = b 2 (9) 

the phase 9 S of b should fall in the interval [0,7r]. It is 
more convenient to instead let 6 S £ [0,2tt], which comes 
at the cost of a local Z-i redundancy. At each site we can 
independently change the sign of both b and N without 
affecting the physical order parameters. As this emergent 
gauge degree of freedom is an artifact of our chosen rep- 
resentation of the operators, any physical operator must 
be a gauge- invariant combination of these operators. In 
particular, inverting Eq. ^ leads to 

b r =s rA /Vv = s r e ie ^ r)/2 , (10) 



where s r depends on a particular gauge, but the combi- 
nation b r s r is gauge invariant, i.e. under a gauge trans- 
formation where b r — > —b r we also have s r — > — s r . 

Any term in the Hamiltonian of b and ./V needs to be 
invariant under such a transformation, in particular the 
kinetic terms must take the form 

-^hopping — tb ^ ^ b r T r r tb r ' ^ ^ N r T r r iN r '. 

{rr') {rr') 

(11) 

Here r , is a Z2 gauge field, which transforms under a 

local Z 2 gauge transformation b r —¥ cr r b r , N r — > a r N r 
as T r y — > a r a r 'T r r ' , with a r G { + 1, —1}- While b itself 
is not gauge- invariant, the integral of its phase 9 S around 
a close loop is, i.e. 

e i2 ^ = {r^ ri r^ r3 ...r^ r )e i ^ (12) 
<j> V0 S)C ■ ds = 2im s , c - (13) 

In the ordered phase (6) ^ 0, single valuedness of b re- 
quires that n s be an integer for any loop. We now con- 
sider a loop enclosing a single defect where 9 C winds by 
2ir, i.e. n c — 1. Then Eq. (JT2J) requires that around this 
loop 

1~r ,ri 7V»2 .rs ■ • • TVjv -V !■ (1^0 

We have seen that single dislocations in e l6c corresponds 
to a Zi gauge flux of tt in the gauge theory. Such an 
excitation (see FIG. fT0| is usually referred to as a vison. 
Therefore the transition between Fermi Liquid and the 
SLM is described as a confinement transition of the Z2 
gauge theory, where visons are gapped in the deconfined 
phase (SLM) while they proliferate in the confined phase 
(FL). At the deconfincmcnt transition, both b and N 
are massive and do not significantly affect the critical 
properties of the gauge theory. However, the phase of 
the gauge theory has important consequences for the b 
and N fields. As long as the visons are gapped, b and N 
are independently well defined. Once the visons condense 
b and N become confined, i.e. they are glued together 
into .^-charge neutral objects like b 2 and bN . 

Let us briefly review some basic properties of Z2 gauge 
theory: A simple Hamiltonian for the Z2 Gauge theory 
on a square lattice is 

// -'IX h Y, II <v (15) 

(ij) □ (ij)ea 

where □ denotes a square plaquette. In addition, the 
system is subject to the gauge constraint 

@i ~ T i,i+x T i,i—x T i,i+y T i,i-y ~ !■ (1^) 
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For h — > (the confined phase of the gauge theory) the 
ground state is 



(ij) (ij) 

In the deconfined phase, J — > 0, it is given by 



|0) 



n 



(17) 

+1)- 
(18) 



Pictorially this state can be represented by coloring in 
gray each bond of the dual lattice which crosses a link 
where tL — — 1 (see FigfTOj). The confined ground state 
is then a superposition of all possible configuration with 
both open and closed gray lines at low energies. The 
vison number is given by the gauge flux n = Y[ n 7-2 thus 
there is a vison at the end of each open line. In the 
deconfined phase (of the Z-i gauge field) the visons are 
gapped, thus there are only closed gray lines. As visons 
(i.e. charge stripe endings) act as sources for these gray 
lines, these lines can be identified with the charge stripes. 



+1 



+1 



T Z = -1 



h > J 



Z2 gauge field deconfined 
Visons gapped 



-EII3- 




J» h 



-1 1- 

.1 1. 



Z'2 gauge field confined 
Visons condensed 



FIG. 10: Left: Typical configuration in the deconfined phase 
of Z2 gauge theory. Electric field lines, = — 1, denoted by 
double black lines in the upper picture, are deconfined. Visons 
- plaquettes with an odd number of ry = — 1 - denoted by a 
spiral in the lower picture, are absent. Visons act as sources 
for the thick gray lines, thus there are only closed strings 
of these. Right: Typical configuration in the confined phase. 
The electric field is confined and visons have proliferated, thus 
there are both open and closed gray lines. In the present 
context the gray lines describe charge stripes. 



C. SLM on tetragonal lattices 

The previous dicussions considered unidirectional 
charge order. On tetragonal lattices there is the addi- 
tional possibility that the order respects the lattice rota- 
tion symmetry, i.e. 



(p(r)) = po 



.Hil 



+ c.c). (19) 



Now there are two possibilities of a Stripe Loop Metal 
respecting lattice rotation symmetries. The SLM can be 
realized by independent fluctuations of the vertical and 
horizontal stripes, i.e. we consider double dislocations in 
8 CjX and 6 CiV separately. Such a state is depicted in FIG. 
[IT] (left). In such a phase there are two kinds of (gapped) 
visons, corresponding to endings of horizontal or vertical 
stripes. 

There is a further possibility for a SLM phase, distinct 
from the previous case, which is illustrated in FIG. [TT] 
(right). This corresponds to allowing, in addition to the 
double dislocations in 9 CX and 9 c _ y separately, a new kind 
of defect where both the phases of both 9 C ^ X and 9 c ^ y wind 
by ±27r. It is easy to see that in the presence of spin 
order, such dislocations are not frustrated. 




SLM (Checkerboard) 



Tetragonal SLM 



FIG. 11: Left: Checkerboard pattern with a double disloca- 
tion in both the vertical and the horizontal stripe. Right: 
There is a different kinds of dislocation, where a horizontal 
stripe turns into a vertical one. Such dislocations also do not 
suffer from spin-frustration, thus there is the possibility of a 
phase where both the double dislocations and this new kind 
of dislocations proliferate. 

In order to describe the transition we begin with two 
equal order parameters (tf) x ) = {i>y) 7^ 0. To include the 
new kind of defect, we write the stripe order parameter 
in a somewhat modified form 



Ipy = 6162,. 



(20) 



Now at a defect where the phase of b\ winds by 27r, both 
9 C X and 9 c . y wind by 2-7T. Similarly, around a single dis- 
location in 62, 9 C X winds by 2-7T while 9 c ^ y winds by — 27r. 
Further, at a single dislocation in both b\ and 62, the 
phase of tp x winds by 4-7T while the phase of ip y does not 
wind at all. Thus this decomposition captures both the 
double dislocations discussed above, as well as the new 
kind of dislocation. Formally this decomposition has only 
a single Zi redundancy, associated with changing the sign 
of both b\ and 62 together. 
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VI. CHARGE STRIPE MELTING 
TRANSITIONS 

A. Orthorhombic crystal 

We begin with the simplest case, i.e. unidirectional 
charge order without spin-order. Our strategy will be to 
initially consider incommensurate stripes (i.e. no pinning 
of the stripe order parameter to the crystalline lattice) in 
the absence of an electronic Fermi surface. We will justify 
this a posteriori by showing that both these couplings are 
irrelevant (in the RG sense). Both in the ordered phase 
and in the SLM the Z2 fluxes (equivalently single disloca- 
tions) are gapped, so we can neglect them when studying 
this transition. To determine the universal properties of 
the phase transition it is convenient to employ a 'soft- 
spin' formulation. The low-energy effective field theory 
is given by the most relevant, symmetry- allowed terms of 
b,N, i.e. 



S[b, N] = / dTd 2 rC b + C N + C 



C b = \Vb\ 2 + -\d T b\ 2 +r b \b\ 2 +u b \b\ 4 



(21) 
(22) 



£ N = \WN\ 2 + \\8 T N\ 2 + r N \N\ 2 + u N \N\ A (23) 



C bN =v\b\ 2 \N\ 2 



(24) 



We note in particular that under a combined time rever- 
sal and spatial inversion (along Q), b is invariant, thus 
no linear derivative terms are allowed. When the spin 
is disordered, N is gapped and may be integrated out. 
What remains is then a 3D XY transition for the b field. 
The critical properties of this model are well knowni^. It 
should be noted that the transition can be formally de- 
scribed as a regular XY transition for the b field, although 
the physical stripe melting transition is in a fundamen- 
tally different universality class. The fundamental field 
b itself is not a physical field, but rather the stripe order 
parameter ip = b 2 . More generally, only operators in- 
variant under the local Z2 transformations are physical, 
thus the operator content of this universality class is very 
different from the usual 3D XY transition. However, the 
critical properties are easily derivable from the ones of 
the 3D XY universality class. The anomalous dimension 
of the composite field^ tp, i.e. r/^ ps 1.49 is dramatically 
different from the anomalous dimension of 6, r\ b ~ 0.04. 
This universality class has been studied in Refs. [33M3^l45l 
where it was dubbed XY* . 

Let us now consider the stability of the XY* stripe 
melting transition to various perturbations. 

1. Pinning to lattice 

The pinning of the period-4 charge stripe to the 
lattice is described by a term 



Spin 



A / drcfr cos(40 c ) 



(25) 



In terms of the phase 9 S of the 6-field this becomes 
Spin = -A / drd 2 r cos(86> s ) (26) 



This is an 8-fold anisotropy on the XY field at the 
2 + 1-D XY fixed point which is well known to be 
strongly irrelevant. Thus despite the commensu- 
rate ordering wavevector the stripe fluctuations de- 
pin from the lattice as the quantum critical point is 
approached. Clearly this pinning term is important 
at long distances in the stripe ordered state, and 
thus represents a dangerously irrelevant perturba- 
tion of the critical fixed point. If the stripe ordering 
occurs at incommensurate wavevector then there is 
no pinning of the stripes either at the critical point 
or in the ordered phase. 

2. Coupling to Fermi surface 

Thus far we have ignored the presence of the con- 
duction electron Fermi surface. It couples to the 
critical stripe fluctuations in several ways which we 
now analyze. The most interesting coupling arises 
in the case where the stripe ordering wavevector 
connects two points of the Fermi surface. In the 
stripe ordered phase this leads to a term 



■k+Q 



Ck 



(27) 



in the conduction electron Hamiltonian which will 
reconstruct the electronic Fermi surface. At the 
critical point or in the stripe melted phase this 
coupling will lead to a damping of the stripe fluc- 
tuations. Integrating out the conduction electrons 
in the stripe melted phase leads to the standard 
Landau-damping term 



\ d / dujd 2 q\uj\\ip(q,uj)\' 



(28) 



This may be viewed as an interaction (between 
the b fields) that is long ranged in imaginary 
time between the stripe fluctuations. The rele- 
vance/irrelevance of this term at the stripe melting 
XY* critical point is readily ascertained by power- 
counting. Under a renormalization group trans- 
formation x — > x' = - , t — > t' = - , we have 



4' -> i/j' = s A tp with A = This implies 



A d — Ad. 



(29) 



As r\ > 1 at the XY* fixed point, the Landau damp- 
ing of the critical stripe fluctuations is irrelevant. 

A direct coupling of the energy density of the stripe 
fluctuations to soft shape fluctuations of the Fermi- 
surface (in the continuum only to the breathing 
mode) is also allowed by symmetry. However, it 
was shown in Ref. that if the correlation length 
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exponent vxy > f , such couplings are irrelevant. 
For the 3D XY model^, v = 0.67155, so this cou- 
pling is also irrelevant. 

For an orthorhombic crystal all other couplings to 
the Fermi surface are even more strongly irrelevant. 
The situation with tetragonal symmetry is more 
complicated and will be analyzed below. 

3. Fermi surface reconstruction 

In the stripe ordered phase gaps will open up at 'hot 
spots' of the original large Fermi surface which are 
connected to each other by the ordering wavevec- 
tor Q. The irrelevance of the coupling of the crit- 
ical stripe fluctuations to the conduction electrons 
suggests that the scale at which the Fermi surface 
reconstructs is parametrically lower than the scale 
at which the stripe order onsets. To explore this we 
now describe the coupled system of stripes and con- 
duction electrons in a different framework that does 
not integrate out the conduction electrons. We re- 
strict attention to conduction electron states near 
the hot spots and linearize their dispersion. The 
resulting action takes the form 



(30) 



S = J drd rC c + C int + C b 
£ c = ^2 Ci {d T - ivi ■ V) Ci 

i 



where Cj denotes a fermion at the ith hotspot, Vi 
is the Fermi velocity at the ith hotspot and tpij is 
a stripe-order parameter with a wave-vector that 
connects the ith and the jth hotspot. Here we gen- 
eralized the fermion action to allow for more than 
one stripe-ordering wavevector, as will be the case 
in the presence of tetragonal symmetry, which we 
will discuss below. Cb is the Euclidean Lagrangian 
of the b held(s) in the 2 + 1-D XY universality class. 

Consider now the RG transformation appropriate 
for the 3D XY fixed point. The action S c = 
J drd 2 rC c is invariant under this RG transforma- 
tion provided that we scale 



Cl,2 



'1,2 = sc l,2 



(31) 



With this scaling the full action is at a 'decoupled' 
fixed point at k = 0. By power counting we see 
that a small n renormalizes as 



KS 



(32) 



As ?7 > 1 k is irrelevant consistent with the anal- 
ysis of the previous section. Now assume we are 
in the ordered phase a distance S from the critical 
point. The energy scale A str i pe of stripe formation 
(i.e the scale at which the critical theory Cb crosses 



over from the critical to the ordered fixed point) 
increases as 



A 



stripe 



1*1- 



(33) 



■^stripe 




/ . g c 
FS reconstruction 



FIG. 12: Stripe ordering and reconstruction of the Fermi- 
surface occur at parametrically different energy scales. The 
dashed lines A str ipc ~ \g — <7c| 0,67 bound the quantum critical 
region and mark the onset of stripe order for g < g c . The 
solid line representing the finite-temperature phase transition 
is parametrically the same as A s tri P e- The dotted line denotes 
the scale where the Fermi-surface reconstructs and is given by 

To determine the value of the gap at the hot spots, 
we consider the self-energy of the electrons. It must 
have the dimensions of energy, i.e. 

£ (w, k, Aslope) = -S (sw, s (j) " 1)/2 k, sA str i p o) . (34) 

We choose s such that sA s tri P o = 1 and scale by a 
factor of K2 Ag tripc /w, obtaining 



£ (U), K, A s tripe 



+773/2 



Astripe' A^T")/ 2 i ' 
\ stripe / 

(35) 



where / is a universal scaling function. On the 
ordered side, close to the hot-spot, the electron self- 
energy has the form (see FIG. IT3|) . 



E(w,fc) 



A 2 

^FS 



(36) 



iu - e k +Q 

where Aps is the size of the gap at the hot spots. 
A FS 

k 

FIG. 13: Diagram for the fermion self-energy in the ordered 
phase. 

By matching this to Eq. (|35|) we get 




A 



(l+'))/2 
Stripe 



(37) 
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Thus as expected the hot spot gap Aps (the scale 
at which the Fermi surface reconstructs) is para- 
metrically different from the scale at which stripe 
ordering occurs. The two energy scales have a ratio 



^stripe 



An-})/2 

Stripe 



(38) 



which goes to zero as the critical point is ap- 
proached. 

Thus the XY* stripe melting critical point survives un- 
modified by either the coupling to the lattice or to the 
electronic Fermi surface. It therefore provides a con- 
crete tractable example of a stripe melting transition in 
a metal. On tuning through this transition the Fermi 
surface undergoes a reconstruction (see Fig. 1141 For a 
more detailed discussion including different parameters 
see Ref. |46| ). This of course happens through the cou- 
pling of the stripe order parameter to the conduction 
electrons as described above in Eqn. [27] 



B. Tetragonal crystal 

We now turn to the situation where in addition to lat- 
tice translations, we also have fourfold rotation symme- 
try. In this case, the average charge density p will be 
modulated in both the x and the y direction, i.e. 



Pr = Pq 



i—x,y 



-2iQi-r 



(39) 



For the translation-symmetry broken state there are now 
two possibilities: Rotational symmetry can be sponta- 
neously broken (ip x ) ^ 0, (tp y ) — (or vice versa), result- 
ing in uni-directional stripes. Alternatively, the charge 
order can preserve rotational symmetry (ip x ) = (ipy) 7^ 0, 
resulting in a "checkerboard" pattern. As we will show 
below, the checkerboard state admits a direct transition 
into the Stripe Loop Metal. In order to reach the Fermi- 
Liquid, both translational and rotational symmetry need 
to be restored. We consider the four possible sequences 
of phase transitions which are shown in Fig. [T5l 

A. Stripe — > Nematic Stripe Loop Metal — > Nematic 
— ► Fermi Liquid 

First we consider a phase transition between a 
unidirectional stripe-ordered state and a nematic 
Fermi-liquid. In this case one of the ipi is zero 
on both sides of the transition. Furthermore, 
rotational symmetry remains broken after trans- 
lation symmetry is restored, so we can describe 
the translation-symmetric state as a condensate of 
pairs of unidirectional stripe dislocations. We thus 
fractionalize the non-zero order parameter as in the 
orthorhombic case and obtain the same action as in 
Eq. 1211 The properties of this phase transition are 
thus identical with the orthorhombic case which we 



discussed in the beginning. After translation sym- 
metry is restored, we still have rotational order and 
the ^-structure remaining in a state which we dub 
Nematic SLM. The ^-structure disappears in a 
second order phase transition^ to a nematic Fermi 
liquid, which is probably first order—. Finally the 
theory of the phase transition between a nematic 
metal and a rotationally symmetric metal 4 ^ has a 
controlled limit 49 , thus the properties of the entire 
sequence of phase transitions are understood. 

B. Stripe -> Nematic SLM^> SLM -> FL 

We again consider the same stripe-melting tran- 
sition into the Nematic SLM, as discussed in the 
previous case. Instead of killing the topological 
structure, the system may first undergo a transition 
where rotational symmetry is restored. This tran- 
sition is unaffected by the presence of the Zi struc- 
ture and is again described by the nematic QCP 
discussed above. The last transition, where the 
topological structure is killed, is slightly different 
than in the case [2U due to tetragonal symmetry. 
This scenario was also studied in Ref. [2f|, and the 
transition is expected to be first order. 

C. Checkerboard-^ SLM -> FL 

We now discuss the case where the ordered phase 
is described by two equal order parameters (ipx) = 
(tpy) ^ 0. To describe a transition out of this state 
we now fractionalize both order parameters 



(40) 



We thus obtain a theory of two order parameters 
coupled to two fluctuating ^-gauge fields. In the 
ordered phase, both horizontal and vertical stripes 
coexist, and defects in the horizontal stripes and 
in the vertical stripes are decoupled. Again our 
aim is to melt the stripes by proliferating pairs of 
dislocations, without closing the gap for single dis- 
locations. In that case we can safely ignore the 
gauge-field and describe the transition by the fol- 
lowing Euclidean action: 



\dM 2 + r\h\ 2 + u\h\ 4 + 2v 



\b x \ 2 \by\ 2 . 



(41) 

We first analyze this action in the absence of any 
coupling to the Fermi-surface and for v = 0. This 
describes two decoupled XY* models which have a 
fixed point where u is of order e = 4 — d, where d is 
the number of spacetime dimensions. The scaling 
dimension of v close to this fixed point is given by 



dim[u] = d — dim[|& x | 2 ] — dim [| 6^1 



v O(2) 



d, (42) 



where vo(2) ~ 0.67155 is the correlation length ex- 
ponent of the AF-model in three dimensions^. We 
thus obtain 



dim[i;] w -0.02 < 0, 



(43) 
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(-27r,27r) 



(27r,27r) (0,tt) 



(7T.7T) (0,71") 



(7T, 7r) 




(tt,0) (0,0) 



for 
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FIG. 14: Reconstruction of the Fermi-surface due to period 4 charge stripes. On the left, the unreconstructed large Fermi 
surface is shown in the extended zone scheme. The dashed lines frame the part of the Brillouin zone that is used to show the 
reconstructed Fermi-surface in the center and on the right. In the center reconstructed Fermi-surface due to uni-directional 
charge stripes is shown in black, while the gray lines indicate the original, folded Fermi surface. On the right, the reconstructed 
Fermi-surface due to checkerboard order is displayed in similar fashion. 




Checkerboard SLM FL 



D 




SLM 



i.e. this perturbation is (weakly) irrelevant and 
the decoupled fixed point is stable. Next, tetrag- 
onal symmetry allows a linear coupling between 
the nematic fluctuations of the order parameter 
ONomatic(r,r) = \b x \ 2 - \b y \ 2 and the ones of the 
Fermi surface 

Pnom. ~ J^(cos k x - cos k y )c\c k , (44) 
k 

where cj. creates an electron with momentum k. 
Unlike the coupling in Eq. ([27j) . this term involves 
fcrmions close to a single point on the Fermi surface 
which leads to enhanced damping |u;| — > |o;|/|g|, i.e. 

ANcraatic J d 2 q J du, J^ I ^Nematic (<?, w) | 2 . (45) 

To determine its relevance, we note that this term 
couples different space-time points, thus its scaling 
dimension is determined by the dimension of | b x \ 2 — 
\b y \ 2 . But at the decoupled fixed point this is given 
by the dimension of \b x \ 2 alone, i.e. 

dim[A N cmatic] = d - 2dim[\b x \ 2 - \b y \ 2 ] = dim[u]. (46) 

Thus the coupling to the Fermi surface is also irrele- 
vant in this case. Therefore the two XY* order pa- 
rameters decouple both from each other and from 
the Fermi-surface at the phase transition and we 
obtain a second order transition from the checker- 
board ordered phase into the SLM phase. 

D. Checkerboard-^ Tetragonal SLM — > FL 

The final sequence we discuss here involves the 
tetragonal SLM introduced in Section IV CI The or- 
der parameters of charge stripes in the x and y 
direction are decomposed as 



FIG. 15: Finite temperature phase diagrams for the three 
sequences of phase transitions which we discuss here. 



4>x = bib 2 
ip y = bib* 2 , 



(47) 
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which comes with a single Z-i redundancy. To de- 
scribe the transition, we thus disorder both ij) x and 
ip y while keeping the gap of the single type of vi- 
son finite. The transformation properties of b\, 62 
under lattice symmetries are shown in TABLE U 



TABLE I: Transformation properties of various fields under 
discrete symmetries (we here assume period-4 charge stripes) 



Symmetry 




i>v 




62 


^-Translation 




ipy 


e i7r/4 6i 


e^ /4 b 2 


^-Translation 


Ipx 


e^y 


e™ /A bi 


e-" /4 b 2 


7r/2 Rotation 


■>Py 


Vx 


62 


bi 


Time Reversal 






bl 


b 2 


x-axis Reflection 






b 2 


h 



On symmetry grounds, this transition is described 
by Eq. (@TJ) and most of the above discussion 
follows through. In particular, the critical stripe 
fluctuations are governed by the same critical ex- 
ponents, and the coupling to a Fermi surface is 
irrelevant. We note that in this representation 
l^zl 2 — IV'yl 2 = 0, thus it appears as if nematic 
fluctuations are strictly zero. However, the opera- 
tor 



Im [ b *i (r + ap)bi (r)] Im [b 2 (r + apt)b 2 (r)] 



(48) 



- ^[i>*(r)tpx(r + afi)]-Re[ip y (r + afi)i> y (r)], 

p=x,y 

(49) 

has the same symmetries as \ip x \ 2 — \"4>y\ 2 - Thus the 
decomposition (|47l) does allow for nematic fluctua- 
tions. 



VII. SPIN STRIPE MELTING TRANSITIONS 

For strong coupling stripe melting transitions it is nat- 
ural to expect that the spin stripe order will melt through 
two phase transitions - first the spin order goes away 
while charge stripe order persists ( i.e translation sym- 
metry remains broken) followed by a second transition 
where the charge stripe also melts. So far we focused ex- 
clusively on the second transition without addressing the 
first. It turns out that this transition (between a spin- 
stripe and a charge-stripe phase) is quite complicated and 
we currently do not have a good theory to describe it. To 
see the source of the difficulty, let us attempt to follow 
the same prescription as above, i.e. initially ignore the 
coupling to the Fermi-surface. Since b is condensed, N 
is now physical (gauge-invariant), thus the transition is 
just the 0(3) transition with tjn ~ 0.04. Now if we add 
the conduction electrons, N can couple directly to their 
spin-density 



k s N ■ ^2 c l + Q&Ck + h.c. 



(50) 
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FIG. 16: Schematic finite temperature phase diagram as a 
function of g (see FIG. [8J . The spin-order vanishes first at g 3 
while the charge-order persists up to g c . 



Following the same argument as below Eq. (|27p we see 
that this is strongly relevant since tjn <C 1 (the marginal 
case is rj — 1), and the fluctuating spin-stripes no longer 
dynamically decouple from the Fermi surface. The re- 
sulting theory was studied by Abanov and Chubuko v 18 i 19 
and by Metlitski and Sachdev^, who showed that it is 
strongly coupled at low energies. There is currently no 
controlled description of this transition. 



A. Stripe multicriticality 

The difficulties with spin-stripe melting can be avoided 
if the spin-melting and the charge-melting transitions 
happen for nearby values of the tuning parameter g (see 
FIG. |HJ), which is likely in the cuprate phase diagram. 
Then, despite the presence of two distinct quantum phase 
transitions the somewhat higher-T physics will be con- 
trolled by a "mother" multicritical point where spin and 
charge stripe order simultaneously melt (see FIG.fT6|). In 
our earlier work-22, we postulated that the temperature 
regime probed in the experiments of Ref. \2]\ is controlled 
by such a multicritical stripe melting fixed point. 

To study this multicritical point where both b and N 
are critical, we again begin by ignoring the presence of the 
Fermi surface and the underlying lattice. The effective 
theory for this transition is the same as in (Eq. 121"]) which 
we reproduce here for convenience: 



S[b, N]= drd 2 r£ b + £ N + £ 



£ b =\\7b\ 2 + - 

|V7V| 2 + 
£ bN = v\b\ 2 \Nf 



-\d T b\ 2 + r b \b\ 2 + u b \b\ 4 



(51) 
(52) 



£ 



■\d T N\ 2 + r N \N\ 2 + u N \N\ 4 (53) 

(54) 



At the multicritical point where both b and N are critical, 
a small coupling v is an irrelevant perturbation^, thus 
we are left with a decoupled 0(3) x 0(2) fixed point. As 
before, we now need to consider the effects of coupling to 
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the Fermi surface and to the lattice. The irrelevance of 
the direct coupling between b and the Fermi surface, as 
well as of the lattice pinning follow by the same argument 
as for pure charge-stripes melting transitions. Like in 
the case of charge stripes, there is a symmetry allowed 
coupling between the energy density of the N fluctuations 
and soft shape-fluctuations of the Fermi surface. Since 
v OCi) > §) this is again irrelevant. Next, the field N is 
not physical (gauge invariant) and thus cannot directly 
couple to the Fermi surface. The gauge invariant spin- 
density that coupled to the spin density of the itinerant 
electrons at the Hot spots is M = bN. In real space and 
(imaginary) time the correlation function of M simply 
factorizes into the correlators of its constituents, i.e. 

(M(r,r)-A?(0,0)) 



l + Vb l + -7s ' 

(r 2 + v 2 t 2 )~ (r 2 + v 2 s t 2 ) — 

In particular, we can read off the anomalous dimension 
of M as rjM = 1 + f]N + Vb- Here 77 at, & are the (small) 
anomalous dimensions of the b and N fields in the regular 
3D XY and 0(3) models, respectively. In the presence 
of a Fermi surface, the fluctuations of the physical spin 
fluctuations M = bN will be Landau damped, like the 
fluctuations of the physical charge fluctuations ip = b 2 in 
Eqn. CM i-e. 



d 2 qduj\uj\\M\ 2 



(56) 



Since t)m > 1, the Landau damping of spin-fluctuations is 
also irrelevant. A further possibility to construct a gauge 
invariant operator as a combination of b and N fields is 

the spin quadrupole operator Q a b = N a Nf, — ^N 2 S a b- It 
couples to the microscopic quadrupole operator of the 
conduction electrons, which is given by 



%(ri,r 2 ) = f(ri - r 2 ) 



a 1 a 2 + cr 1 cr 2 



(57) 



where ai is the spin-density at rt and f(r) is a non- 
universal function obeying lattice symmetries. For slow 
fluctuations of Qrg all the electrons need to be on the 
Fermi-surface which imposes strong phase-space con- 
straints familiar from Fermi-liquid theory. It is thus con- 
venient to decompose the vertex into its angular harmon- 
ics in the particle-hole and Cooper channel, to get 



(58) 



Jo ft m 

/ dujY,Y. V ^ C '^ LJ ) C ~ra{K-k i n-L0) 1 

Jo „, 



k m 



where Tl m and C rn are the particle-hole and particle- 
particle propagators, respectively, with angular momen- 
tum m. Now it is straightforward to evaluate XQ term by 



term and thus establish the irrelevance of Landau damp- 
ing. In particular in the m = channel we have, for a 
Fermi-liquid 



C^(uj,k) - 6(o; 2 -v 2 F k 2 ) 



y/(v F k¥ 



(59) 
(60) 



and it is then clear that x'q ~ ^ 3 (with a cut-off depen- 
dent logarithm in the particle-hole channel) and hence 
damping is strongly irrelevant. 

We have thus verified that the decoupled multicriti- 
cal point is indeed stable both against coupling to the 
Fermi surface, and against pinning to the crystalline lat- 
tice. As we pointed out, this is the opposite of what 
one would expect for a conventional stripe melting tran- 
sition, where the stripe fluctuations become strongly cou- 
pled to the Fermi surface at low energies. This immunity 
of the critical stripe fluctuations to the properties of the 
Fermi-surface, whether it is gapless or has a single par- 
ticle (pseudo-)gap is consistent with the experiments of 
Ref. [23. More quantitatively, in the theory introduced 
above, stripe correlations at finite temperature obey oj/T 
scaling with a dynamical critical exponent z = 1, as ob- 
served experimentally. Further, the height of the peak 
in the dynamical spin-susceptibility scales with temper- 
ature (see Sec. IX Aft as lim — ,^"' T ' ) ~ , where 



r\M = 1 + Vn + Vb ~ 1-08. This again agrees remarkably 
well with the results of Ref. 



VIII. SINGLE PARTICLE GREEN-FUNCTION 
AND FERMI SURFACE RECONSTRUCTION 



We now turn to calculate the Green-function of the 
electrons at the phase transition. Since the coupling to 
the fluctuating stripe-order is irrelevant, it is sufficient 
to consider the leading order contribution to the electron 
self-energy in perturbation theory. The same self-energy 
has already been evaluated explicitly in a different con- 
text in Ref. l33l . We here reproduce the universal behavior 
from scaling. Electrons on points of the Fermi-surface 
which are connected by Q (hot spots) will be affected 
most strongly, so we focus on those. In the vicinity of 



k k — q k 



FIG. 17: Diagram for the leading contribution to the Fermion 
self-energy. 



the critical point, the electron self-energy obeys (see Eq. 
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v ,, v ,'dk z A stripe k 
S (w, dfc, k, A stri p C ) = ^5 , , 



(61) 



where <5fc is the momentum-deviation from the hot spot 
and g is a universal scaling function. The leading order 
contribution (see FIG. IT?]) appears at order k 2 , so we get 



£(w,(5fc = 0,K^0,A stri pc = 0) 



(62) 



Since t] > 1 the Fermi-surface as well as the Lan- 
dau quasiparticles remain sharply defined at this critical 
point, even at the hot-spots. 

We can further use the scaling form of the self-energy 
to discuss the reconstruction of the Fermi surface. For 
sufficiently large doping g > g c the stripes are melted 
and we have a large Fermi surface. At lower doping 
the stripe fluctuations, and consequently the scattering 
of Fermions increases. It is clear that the scattering is 
non-uniform on the Fermi-surface, i.e. it decreases with 
distance to the hot-spot. At even lower doping, spectral 
weight from the original, large Fermi-surface starts being 
transferred to a "shadow" Fermi surface, i.e. the surface 
of k points that is connected to the large Fermi-surface by 
Q. Deep in the stripe ordered phase, the shadow Fermi 
surface becomes part of the new, reconstructed Fermi 
surface. In the vicinity of the critical point, the spectral 
weight obeys universal scaling (with non-universal angle- 
dependent prefactors) which can be easily calculated. 

7 

^ Shadowband 
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FIG. 18: The spectral weight on the shadowband vanishes as 
a power law close to the critical point. 

We compute the self-energy for electrons (FIG. H7|) 
close to a generic point Kq on this "shadow" Fermi- 
surface, safely away from the hot spots. It is again given 
by Eq. (jrTTj) . since the intermediate electron lies on the 



Fermi surface. For cj/vf, Sk, A^ripe/Vp "C K the de- 
nominator in the imaginary part of the Green function 

lmG(w, K + 5k, A str ipo) = 



(w - e Ka+Sk - E') 2 + (£") 2 
(63) 



can be replaced by e 2 K and the spectral function inherits 
the scaling form of E , i.e. 



A(u,K + Sk) ~ -j— G( 



Sk z A. 



stripe v 



^K 



U) CO 



where G is a new scaling function. 

This implies^ that the quasiparticle weight Z. 
obeys a power law (see FIG. [T5)l 

Shadow ~ ^-\g - g c \ u{z+n) ■ 

6 K„ 



(64) 



shadow 



(65) 



Here 77 is the anomalous dimension of either charge of spin 
fluctuations, depending on which ordering is responsible 
for the reconstruction. 



IX. PAIRING BY STRIPE FLUCTUATIONS 

We now turn to the possibility of superconductivity 
due to stripe fluctuations. We first focus on the stripe 
melted phase and the approach to the critical point. 
Since the coupling between the fermions on the Fermi- 
surface and the fluctuating stripes is irrelevant, it is safe 
to integrate out the stripe-fluctuations which gives rise 
to a 4-fcrmion interaction 



V(co,k) = 



Q-k\ 2 ~^ 



■WQQ-k\ZM), 



(66) 



where PL is a scaling function with. 

lim W(x, 0) = x 2 - Vs lim W(x, 0) = 1 (67) 

x— >0 x— >oc 

A sample function with these properties is 

x 2-r, s 



W(x,G) 



c 2- Vb + I 



(68) 



This interaction can then be treated in mean-field theory. 
Since the interaction is repulsive, it can only contribute 
to a superconducting order parameter with opposite sign 
on hot-spots which are connected by Q s (see FIG. [H?)) . 
For Q s close to (it, n) and a typical (large) Fermi-surface 
of the cuprates, the lowest angular momentum channel 
that satisfies this is^ d x 2_ y 2. To determine the mean- 
field gap at zero temperature, at weak coupling, we can 
neglect the frequency dependence of the interaction, i.e. 



exp 



1 



N{E F )u d \ 
dkdk'V{k - k')d k d k .8{e k 



(69) 
(70) 
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FIG. 19: (color online) Left: A typical large Fermi-surface of 
the cuprates is shown in the first Brillouin zone, where the 
shaded regions denote the sign of a d x 2_ y 2 order parameter. 
The solid dots denote a pair of electrons, which resides at 
momenta where the sign of the order parameter is positive. 
After scattering on a fluctuating stripe, the scattered Cooper 
pair, denoted by the tips of the arrows, find itself at momenta 
where the sign of the order parameter has changed. Since the 
interaction itself is repulsive, the effective electron-electron 
interaction is attractive in the d x 2_ y 2 channel . Right: Phase 
diagram including superconductivity induced by stripe fluc- 
tuations. The transition temperature T c is highest when spin 
fluctuations are critical. 



where dk is a function with d x i_ y i symmetry. The dom- 
inant contributions to this integral come from momenta 
close to the hot-spots where dkd/~' ~ — 1. From these 
momenta we get 



X. EXPERIMENTS 

In this section we discuss our results for the quan- 
tum critical stripe fluctuations in the context of available 
data. 



A. Scaling of quantum critical spin fluctuations 

We first briefly review the standard scaling assump- 
tions for the full q and ui dependent dynamical spin sus- 
ceptibility associated with quantum critical spin fluctua- 
tions. The imaginary part can be expressed in terms of 
a universal scaling function F as 



x"(q,u;T,6) 



co 



-F 



u d ^,K)^ K 2 C Vs U{iK), 



(71) 



\q-Q s \ 2 ^ \ Cl \q-Q s \^T' E S" 

(75) 

As before, Q s is the incommensurate peak wave vector, 
T is temperature and 5 — \g — g c \ is the distance from the 
T = quantum critical point, z is the dynamical critical 
exponent, v is the correlation length exponent, and rj s 
is the anomalous dimension of the spin order parameter. 
Co, ci, Eq are non- universal numbers. 

First, at to, T 3> Eq5 uz the system is in the "quan- 
tum critical regime" and does not know that the ground 
state is not exactly at the quantum critical point. In this 
regime the width of the peak in g-space will satisfy 



where K is a momentum cut-off and the scaling function 
U must obeys 



lim U (x) = x^ 3 

x^oo 



lim U(x) = x 2 . (72) 

x— >0 



Next we want to determine the transition temperature 
T c . It is clear that for temperatures T 3> the fre- 
quency dependence of the interaction can no longer be 
dropped. For a qualitative result it is sufficient to im- 
plement this approximately by writing the temperature 
dependent effective interaction as 



u d {£, K- T) a Ud (^e+T- 2 , K) 



(73) 



For weak coupling, the self-consistent equation for T c 
then becomes 



1 A _ 1 

° S T c u d {T c )N{E F ) ■ 



(74) 



where A is a UV cutoff. 

In the stripe ordered phase, the region around the hot- 
spots becomes gapped. This will suppress the ability 
of the residual stripe flutuations to cause pairing. The 
resultant phase-diagram is shown schematically in FIG. 
[PHI As expected we find that pairing is strongest close 
to the critical point, where it is mediated by gapless spin 
fluctuations. 



k(u,T) 



ci 



T 



(76) 



where K is a universal scaling function. Restricting to 
u> = then we get the well known result 



k{uj = 0,T> E„S UZ ) ~T». 



(77) 



Second, exactly at the peak q = Q, we have (still in 
the "quantum critical" regime) 



(78) 



with ao a non-universal number (related to Co,Ci) and 
X a universal scaling function. At low frequency at a 
non-zero T, \" is linear in u) so that the scaling function 
X{x) ~ x for small x. This then gives 



lim 



Xp(w,T» E^) 



«,() 



T{2- Vs + z)/z 



(79) 



B. Scaling in the LSCO experiment of Ref. l27l 

We will now use the scaling properties reviewed above 
to carefully discuss the experiment in Ref. H3- There, 
Aeppli et. al. used neutron scattering to find the spin 
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structure factor S(Q,uj; T) which is related to the imag- 
inary part of the susceptibility via 

S(Q,u);T) = — — > . 80) 

1 — e t ui 

In their data, they fit the temperature dependence of the 
peak height to a power-law (Fig. 3B in Ref. H3 ) 



Hm XP.Fit(^>^) „ T - 1 .94 ) 
w— >0 UJ 



from which it follows (|79l) that 

2 - T) s + z 



1.94. 



(81) 



(82) 



To determine z, the width of the peak n is plotted as 
a function of \fus 2 + T 2 over a range of frequencies and 
temperatures (Fig. 4 in Ref. 27 ). A linear relationship 



k{uj,T) ~ yjuj 2 + T 2 



(83) 



is found, which implies (|76|) that z = 1, which in turn 
determines the anomalous dimension r\ a s» 1.06. 

Some caution is however required with the quantum 
critical interpretation of the data. To emphasize this we 
note that Ref. |27J also plotted Yp(ui, T) as a function of 
k(uj, T) (Inset of Fig. 4 in Ref. (27| ) and from a fit to the 
data found that 



lim 



X ^ Fit (u,T^>Eo6" 



k~~ 3 (lo =0,T). (84) 



If quantum critical scaling applies, then from Eq. (|77l79p 
it follows that 



lim 



x';^,t>e ^) 



UJ 



K n.-3( u = 0,T). (85) 



Thus the last fit implies r; s = 0, in apparent contradiction 
to the conclusion rj s — 1.06, reached from the previous 
fits. 

To address this question, it is necessary to analyze 
k(lo,T), which relates the two apparently contradictory 
fits. k(lu,T) saturates at low T,uj. One possible cause is 
of course that the system is not exactly at the quantum 
critical point. A second possibility is quenched disorder. 
As the stripe order parameter breaks translation symme- 
try, non-magnetic disorder will couple to it roughly as a 
"random field" . This will always lead to a saturation of 
k - even when there is true stripe ordering. 

One resolution of this apparent contradiction is that 
the fit leading to (jMl) includes data from the region where 
k has saturated. If the saturation is mainly due to disor- 
der, those data points (small values of k) should not be 
used in scaling plots. The remaining data points can be 
fit with 



lim 



X 'k Fit (u,T ^ E 6»*) 



k- 2 (w = 0,T), (86) 



in agreement with r/ s w 1.06. Of course, once several data 
points are excluded, the range of the remaining data is 
rather small and consequently the error in extracting the 
exponent is large. 

In a more recent paper—, the same neutron data is 
used to compute 



T c ff = T / lim 



X"(<7,^,T) 



(87) 



If the q integral is restricted to the scaling region around 
Q, scaling applies and the contribution of those wave- 
vectors to r e ff is given by 



7cff,scal — T 



/ lim 

JqnQ 



X"(9,w,T) 



Thus, since rj s is close to one, T e ff iSC ai must have a nearly 
linear T dependence. In the data of Ref. [H, r e g is 
roughly T-independent. If scaling is indeed satisfied by 
the neutron data near Q then it must be that the q- 
integral is dominated by non-scaling contributions away 
from the peak. 

These caveats about the evidence for critical scaling in 
the existing data strongly emphasizes the need for further 
experimental efforts to study the singular stripe fluctua- 
tions in modern samples of cuprates. Little is also known 
about the doping and magnetic field dependence of these 
singular fluctuations. In light of the suggested crucial 
rolei^ of these fluctuations for the physics of the strange 
metal and the theory discussed in this paper we hope 
that such experimental efforts will be forthcoming in the 
near future. 



C. Other scaling in cuprate neutron experiments 

Scaling of the spin fluctuation spectrum has been re- 
ported in a number of cuprates over the years^~— . In 
an early paper, Keimer et al" studied very lightly doped 
LSCO and found evidence of scaling of the local (i.e q- 
intcgrated) dynamic susceptibility as a function of the ra- 
tio uj/T. Subsequent experiments on highly underdoped 
YBCO near the edge of the superconducting boundary 
also saw similar uj/T scaling 55 . This scaling in very un- 
derdoped cuprate samples presumably has little to do 
with the stripe quantum criticality discussed in this pa- 
per. 

More pertinent possibly is the uj/T scaling seen^ in 
the normal state of underdoped YBCO6.5. This is in 
the pseudogap regime and the spin fluctuation spectrum 
has incommensurate peaks ( "dynamic stripes" ) . Detailed 
studies of the q and T dependence of the scale invariant 
part of the scattering have not been reported in the liter- 
ature. In the context of this paper it is particularly inter- 
esting to ask if this scaling is associated with the 'mother' 
multicritical regime (see Fig. I16[) at intermediate tem- 
peratures. If so then the detailed q, ui and T dependence 
of x" will satisfy the scaling form discussed above with 
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the exponent r/ s 1.06, z = 1. Upon cooling supercon- 
ductivity develops and the measured spin fluctuations 
change character. It is interesting to ask what happens if 
the superconductivity is suppressed in a field. At similar 
doping levels in YBCO recent NMR experiments^ show 
that there is charge stripe order without any accompany- 
ing spin stripe order. In terms of Fig. [Tj)]this places the 
system between g c and g s . In this regime in our theory 
there are dynamic scale invariant stripe fluctuations at 
intermediate temperatures. Upon cooling the spin cor- 
relation length saturates to a finite value but there is a 
charge ordering transition. Associated with this there is 
a reconstruction of the Fermi surface. 



D. Dynamic susceptibility 

For comparison with future experiments, we now com- 
pute the dynamic spin-susceptibility at the multicritical 
point, i.e. the Fourier-transform of Eq. (|55l) : 



vi\K\ 



-v 2 2 {k 



(89) 
K) 2 ) 



Two trivial limits of this function are given by 



e(n 2 - vlK 2 ) 



(n 2 - v 2 k 2 ) 



1-771-772 



2 -772 
2 

(90) 
(91) 



For finite K and V\ ^ i>2 one can see that the most 
singular contributions to x'li{n,K) occur when all the 
momentum and frequency are carried by either b or N, 
which is possible only for VL = v\K or Q = v^K. For ex- 
ample if Q — v%K then the integrand diverges at w, k — > 
as ui^ +r,1 ~ 3 . After integration a power-law singularity 



xm^k + en, k) - xlAv2K, k) ~ sn~ 



(92) 



remains. A numerical plot of this function is shown in 
FIG. 1201 In any actual experiment, xm will be modified 
from this by finite temperature effects. In particular, 
all sharp features (onset and peaks) are expected to be 
washed out. We expect that what remains as a qualita- 
tive feature is a dispersive peak which is relatively narrow 
as a function of K but has a long tail as function of O. 




FIG. 20: The imaginary part of the spin-susceptibility, plotted 
for three fixed values of Q, as a function of K (MDC). The 
value of x"(^j-f) is given by the height of the blue curve 
above the (dotted) base line. For any given value Q it is non- 
zero only over a finite interval. For fixed K and as a function 
of fl (EDC) it has a similarly sharp onset at v±K, but does 
not drop back to zero at higher frequencies. 

Sec. IX Al In particular, right at the critical point, the 
peak amplitude at low frequencies obeys 

Sp,Chargc(T) ~ (93) 

where rj c is the anomalous exponent of charge fluctua- 
tions, which is rjc 1.49 in our theory. Measuring this 
exponent thus constitutes an independent test of our the- 
ory. 



F. Tunneling density of states near pinned stripes 

STM has been a highly successful tool in investiga- 
tion of striped phases in the cuprates. Recently, it has 
been used to image individual defects in the stripe order 
of BSCCO 57 . Within our theory as the stripe melting 
transition is approached the core energy of doubled dis- 
locations decreases to zero while that of single ones stays 
non-zero. It is possible therefore that upon approaching 
optimal doping the density of close dislocation pairs in- 
creases at the expense of isolated dislocations. It will be 
interesting to explore this in STM experiments. 

STM experiments can also potentially be used to deter- 
mine the anomalous dimension of the charge stripe fluc- 
tuations i] c . In the stripe- melted phase, but close to the 
phase transition, there will always be disorder which lo- 
cally pinns fluctuating stripes. Close to such an impurity, 
there is then local stripe-order with an amplitude which 
decays with distance from the impurity. By locally mea- 
suring the tunneling density of states, these oscillations, 
as well as their decay, can be detected experimentally. 



E. Scaling of quantum critical charge fluctuations 



In addition to measuring critical spin fluctuations via 
neutron scattering, it may be possible to detect criti- 
cal charge fluctuations, e.g. via X-ray experiments. At 
a quantum critical point where charge order melts, the 
scaling properties of the dynamical charge structure fac- 
tor Scharge(Q, T) follow from the same analysis as in 



FIG. 21: The lowest order contribution of pinned stripes to 
the electronic tunneling density of states. The box denotes a 
time-independent disorder potential 
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FIG. 22: Density of states for tunneling in electrons at zero 
bias as a function of distance from an impurity. The impurity 
locally pins stripe fluctuations which lead to an oscillatory 
charge density close to the impurity. Conduction electrons 
scatter on this pinned stripe fluctuations which leads to an 
oscillatory tunneling density of states. 

The lowest order oscillatory contribution to the tun- 
neling density of states is shown in FIG. [5TJ The con- 
tribution due a delta-function disorder potential at the 
origin is 

N u (x)~ / G x - v (uj)cos(Q ■ y) D y (n = 0)G y - x (uj) 
Jy 

(94) 

We note that the integral is dominated by electrons 
close points on the contour of energy /x + U) which are 
connected by Q. Around these points we can expand the 
dispersion to linear order, i.e. under the integral we have 
Gui (x) ~ - . The propagator of charge-stripe fluctuations 
at zero frequency is D(x) ~ — ^ so we get 

N u=0 (x) ~ Iccr^ cos (Q + (95) 

where $ is some constant. The spatial decay in the am- 
plitude of the quasiparticle oscillations is directly given 
by the anomalous exponent rj c of the charge stripes (see 
FIG. HI]). We note that unlike quasiparticle interference, 
these oscillations are non-dispersive, i.e. the wave-vector 
Q is set by the stripe order, and does not change as a 
function of bias uj. 



XI. CONCLUSIONS 

In this paper we have studied several examples of con- 
tinuous quantum melting phase transitions of stripe or- 
der motivated by the cuprates. The stripe melting was 
assumed to be driven by the proliferation of doubled dis- 
locations. We developed a simple and intuitive picture of 
the resulting stripe melted phase: the fluctuating stripes 
form closed loop configurations of arbitrary size. Cut- 
ting a loop open to produce two end points costs finite 
energy. We therefore dubbed this phase a "Stripe Loop 
Metal" . A Stripe Loop Metal phase has sharply defined 



Landau quasiparticles at a Fermi surface that satisfies 
the usual Luttinger theorem. It however differs from the 
conventional Fermi liquid in some very subtle ways. The 
spectrum of excitations of the stripe order parameter is 
fractionalized. Associated with this there is a gapped 
Z2 topological defect that is the end point of an open 
stripe(equivalently a relic of the single stripe dislocation 
of the ordered phase) . These differences are sufficiently 
delicate that even if such a stripe melted state exists in 
a material, it cannot easily be distinguished from a con- 
ventional Landau Fermi liquid. 

We showed that this kind of stripe melting transition 
is strongly coupled yet tractable. The critical stripe fluc- 
tuations dynamically decouple from the Fermi surface. 
However, in the ordered phase, the stripes recouple to 
the Fermi surface and reconstruct it. Thus our results 
provide concrete examples of tractable strongly coupled 
quantum phase transitions associated with stripe melt- 
ing and the associated reconstruction of the Fermi sur- 
face. We are not aware of any other controlled theories of 
quantum criticality related to Fermi surface reconstruc- 
tion by translation symmetry breaking order in two space 
dimensions. 

In the examples studied in this paper the dangerous ir- 
relevance of the coupling of the stripe fluctuations to the 
Fermi surface implies that the latter reconstructs at a 
scale that is parametrically lower than the scale at which 
the stripe ordering occurs. Furthermore the stripe fluc- 
tuations do not destroy the Landau quasiparticle at any 
point of the Fermi surface at the quantum critical point. 

Before concluding we discuss these results in the con- 
text of several recent experiments and theoretical sugges- 
tions on the cuprates. As discussed in the introduction 
there is growing evidence for the ubiquity of stripe order 
in the underdoped cuprates. Stripe order has also been 
invoked to reconstruct a large Fermi surface to obtain 
pockets which may explain quantum oscillation phenom- 
ena in a magnetic field at low-T, and Hall effect and 
other measurements. In contrast there is very little evi- 
dence for stripe ordering in the overdoped side. A natu- 
ral suggestion is that a quantum critical point associated 
with melting of stripe order and the related Fermi sur- 
face change underlies some of the strange normal state 
properties around optimal doping. 

To apply the results of this paper to the cuprates we 
first need to postulate that the overdoped metal is a 
Stripe Loop Metal. There is very little experimental evi- 
dence against such a proposal, and so it is worthwhile to 
examine our results further. As proposed in our earlier 
work>22 and discussed further in this paper, the 'mother' 
multicritical fixed point where the charge and spin stripe 
orders simultaneously quantum melt into an overdoped 
Stripe Loop Metal might potentially control the physics 
of the stripe fluctuations around optimal doping in the 
strange metal regime. This proposal finds some support 
in well known early experiments of Aeppli et al^I on near- 
optimal La2- a; Sr 2 ;Cu04 providing an explanation of both 
the scaling of the peak width and peak height as a func- 
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tion of temperature/frequency. We outlined a number 
of other predictions from the proposed theory for future 
experiments. 

A disappointing aspect of our results is that the crit- 
ical stripe fluctuations associated with the transition to 
the Stripe Loop Metal are not capable of producing most 
of the observed non- Fermi liquid physics of the strange 
metal regime. In particular the electron quasiparticle re- 
mains well defined even at the stripe melting quantum 
critical point. Given the success of our theory in describ- 
ing the observed stripe fluctuations in neutron data one 
possibility is simply that the physics driving the strange 
metal non-fermi liquid is not stripe quantum criticality 
(contrary to the proposal of Ref. 0). We suggest that 
our theory correctly describes the stripe sector which de- 
couples dynamically from the electronic fluctuations of 
a non-fermi liquid metal at the stripe quantum critical 
point. The explanation of the non- Fermi liquid physics 
itself should then be sought in some other mechanism. 

We thank G. Aeppli, J. C. Davis, Eduardo Fradkin, 
Eun-Ah Kim, Steve Kivelson, Michael Lawler, Young Lee 
and Stephen Hayden for useful discussions. TS was sup- 
ported by NSF Grant DMR-1005434. 



Appendix A: Landau-damping in some non-Fermi 
liquid metals 

Consider a metal with a sharp Fermi surface, possibly 
in a Non-Fermi Liquid phase. We take a scaling form for 
the electronic spectral function 



A(8k,w) = 



1 



y./ Z 



(Al) 



where 8k = (k — kp) ■ kp. This is appropriate for the 
Fermi-Liquid (a = Zf = 1) but also for the marginal 
Fermi-Liquid and certain non-Fermi Liquids, such as a 
metal at a Pomeranchuk transition (a — 1, zj = 3/2). 
Using this to perturbatively calculate the self-energy of 
the fluctuations of the stripe order parameter (FIG. [3]), 
we find 

U"(u,Q) = J dfl J d 2 kA(k,u + n)A(k,Lo) (A2) 

^uj 1+2lT f. (A3) 

At low energies, for a > 2 2 Z/ , this dominates over the 
bare w 2 term and consequently z is renormalized away 
from one. 

This remains true even for some more exotic models 
of the strange metal, such as the uniform Resonating 
Valence Bond (RVB) state studied, e.g., in Ref. |58|. In 
this model there is spin-charge separation, with the spin 
carried by charge-neutral fermions (spinons) which form 
a Fermi surface, and the charge carried by spinless bosons 
(holons) which form an incoherent Boltzmann gas. The 
electron Green function is given by the convolution of 



a spinon and a holon Green function. Nevertheless, at 
finite wave-vectors the Landau-damping again takes the 
form of Eq. (|A3j) . The spin response of this state is 
determined solely by the spinons (which form a Fermi 
surface), thus the self-energy of spin stripe fluctuations 
has the form of Eq. (|A3p . The charge response is given 
by the Ioffe-Larkin rule 



n 



n 



electron 



spinon-H-holon 



(A4) 



^spinon 



where the holon-polarizability is finite and non-singular 
at finite wave-vectors and LT" pinon is of the form of 
Eq. (|A3|) . Expanding n" lcctron for small ui then yields 



electron 



Appendix B: RG-equations for the Tetragonal model 

Consider the theory described by Eq. (14"TT) . i.e. 

J T T J T r 

z — x,y ' ' ' ' ' ' 

(Bl) 

where we wrote ipi = (Re[6j],Im[6j]). More generally 
we can consider this theory for TV-component vectors ipi . 
The RG-equations to leading order in e = 4— d are readily 
obtained 



dl 
d_ 

Jl 



u = eu-(8 + N)u 2 I - Nv 2 I a (B2) 
v = ev-{4 + 2N)uvI + 4:V 2 Io, (B3) 



where Iq depends on the cut-off scheme. The fixed 
points (FPs) are 

A : 
B : 



C 
D 



u , v 



=(0,0) 



; + 7V 
N 



,0)e// 

N-4 



16 + 27V 2 ' 16 + 2N 2 
1 1 



e/Io 



8 + 27V' 8 + 2N 



e/ io- 



(B5) 
(B6) 
(B7) 



A standard stability analysis finds that the 0(27V)-fixed 
point D is stable for N < 2, C is stable for 2 < N < 4 
and the decoupled FP B is stable for N > 4. However it 
is known from higher-order expansions and Monte-Carlo 
studies that for N > 1 the decoupled FP is the only stable 
fixed point in d = 3. So the leading order e-expansion 
gives the incorrect result. Similarly, the leading order 
expansion ml/N finds the decoupled FP to be unstable 
(this can be seen from Eq. (|4"2"|) , where to leading order in 
1/7V, v = 1/2). A higher order expansion in the presence 
of the nematic coupling (j45|) . e.g. along the lines of Ref. 
[25l . is rather involved, and we do not pursue this direction 
here. 
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Appendix C: Dual description of the deconfined 
transition 

First consider the orthorhombic case. The dual de- 
scription of the usual AY-transition is given by a vortex 
field ^(r) = e lv ^ which is minimally coupled to a fluc- 
tuating U(l) gauge field, i.e. the hopping of vortices is 
described by 

H v = —t v cxp{iip(r) — iip(r + e) — 2nia e (r)} + h.c, 

r,e 

(CI) 

where the number of the original AY-bosons is related 
to flux of the gauge field via 

n r = cfi(r + x) — cfi(r) - a x (r + y) + a x (r) = V x a. 

(C2) 

In the following we will use V to describe lattice deriva- 
tives to lighten the notation. When the vortices are con- 
densed (&) 7^ 0, single valuedness of * requires that 

2ir i ■ a = 27r V x a = 2nN, (C3) 

t&n 

i.e. the boson number is quantized in integers. For the 
XY* transition, we only allow doubled vortices $(r) = 
e 2i<p(r) _ Their hopping is described by 

H* = -t v exp{i2Vip{r) - Ama{r)} + h.c, (C4) 

r 

and when ($) 7^ 0, single valuedness of $ implies 

4ttV x a = 2ttN, (C5) 

i.e. the boson number is quantized in half-integers, in 
agreement with our fractionalizion prescription 

ip = b 2 . (C6) 



Further, the single-vortex is now an Ising variable, i.e. 
$ = ±y/ ($) and can be identified with the Z2 gauge 
field. 

Now we turn to the tetragonal case, where we have two 
vortex fields ^i(r) = e l f^ r \ each coupled to their own 
gauge-field at(r). We only allow hopping of pairs ^1^2 
or *i*2i i- e - 



Hl ct = -t v Yl exp{iV^(r)+aiV^ 2 (r) (C7) 
- 2m{a\{r) + aa 2 {r))} + h.c. 

If (*i* 2 ) = (*i*2> the fluxes are quantized as 



Vxo 1 ±Vxa 2 =iV ± , (C8) 



or alternatively the boson numbers as 

ni = i(JV + +JV_) (C9) 

n 2 = i(iV+-iV_), (CIO) 

so in terms of the charges Qi,Q2 of the original XY- 
fields, the elementary excitations in the phase where the 
paired vortices have condensed carry charge (Qi,Q 2 ) — 
(1/2, ±1/2), which is reflected in the fractionalization 
prescription (Eq. ([2010 



ipx = bib 2 

'4>y = hb* 2 . (Cll) 
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